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Abstract. Recently, we have derived a generalization of the known canonical fluc- 
tuation relation between the heat capacity and the energy fluctuations C — (3 2 (SU 2 ), 
which is able to describe the existence of macrostates with negative heat capacities 
C < 0. In this work, we extend our previous analysis to the generalization of some 
other equilibrium fluctuation relations such as \t = ft (<5M 2 ) and Kt = P {^V 2 ) /V, 
in order to consider the existence of macrostates with anomalous values in response 
functions such as the isothermal susceptibility in a paramagnetic system \t = BM/dH 
or the isothermal compressibility of a fluid system Kt = —V^ 1 (dV/dp) T . Finally, we 
discuss some simple application examples and perform a numerical testing of some con- 
sequences of the present generalization by studying the anomalous diamagnetic states 
of the 2D Ising model through Monte Carlo simulations. 
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1. Introduction 

According to the well-known fluctuation relation: 

C = k B (3 2 (5U 2 ) (1) 

between the heat capacity C and the energy fluctuations (5U 2 ), the heat capacity should 
be nonnegative. However, such a conclusion is only an illusion. Since the first theoretical 
demonstration about the existence of macrostates with negative heat capacities C < by 
Lyndel-Bell in the astrophysical context PQ, this anomaly has been observed in several 
systems [21 [31 IU El El [7J [8]. The fluctuation relation directly follows from the 
consideration of the Gibbs' canonical ensemble: 

dpc (U\(3) = ^ exp (-PV) n (U) dU, (2) 

which accounts for the equilibrium thermodynamic properties of a system in thermal 
contact with a heat bath at constant temperature T when other thermodynamical 
variables like the system volume V or a magnetic field H are kept fixed, where 
(3 = 1/kgT. As elsewhere shown [H [9], macrostates with C < can exist within the 
microcanonical description of a given system, but these states are thermodynamically 
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unstable under the external influence imposed to the system within the Gibbs' canonical 
ensemble. 

A suitable extension of this standard result has been recently derived [TO] : 



which considers a system-surroundings equilibrium situation where the inverse 
temperature [5 W of a given thermostat exhibits correlated fluctuations with the total 
energy U of the system under study as a consequence of the underlying thermodynamic 
interaction. This last generalization differs from the canonical equilibrium situation in 
the fact that the internal thermodynamical state of the thermostat can be affected by the 
presence of the system under study, constituting in this way a more general framework 
than the usual canonical (<5/3 w = 0) and microcanonical ensembles (5U = 0). 

The new fluctuation relation (j3J) defines a criterion capable of detecting the presence 
of a regime with C < in the microcanonical caloric curve 1/T(U) = OS (U) /dU of 
the system through the correlated fluctuations of the inverse temperature /3 W and the 
energy U of the system itself. In fact, it asserts that macrostates with C < are 
thermodynamically stable provided that the influence of the thermostat obeyed the 
inequality (5/3 u 5U) > 1. Consequently, any attempt to impose the canonical conditions 
8fi u — > leads to very large energy fluctuations 5U — > oo that induce the thermodynamic 
instability of such anomalous macrostates. 

As already discussed in a previous work [TT], the fluctuation-dissipation relation 
03]) has interesting connections with different questions within statistical mechanics, 
such as the extension of canonical Monte Carlo methods to allow the study macrostates 
with negative heat capacities and to avoid the super-critical slowing down of first-order 
phase transitions [12], the justification of a complementary relation between energy 
and temperature p31 HU [151 EE CEO US CHI EQl EJJ, the development of a geometric 
formulation for fluctuation theory [221 [23] , as well as the Mandelbrot's derivation 
of statistical mechanics from inference theory [16] . However, Eq.([3]) is applicable to 
those equilibrium situations where there is only involved the conjugated pair energy- 
temperature. Thus, this result merely constitutes a first step towards the development 
of an extension of equilibrium fluctuation-dissipation relations compatible with the 
existence of anomalous response functions [271 1211 [29], [301 ED [32] . 

In this paper, we shall extend our previous results to an equilibrium situation with 
several control parameters, that is, we shall develop a special approach of fluctuation 
theory [24] in order to arrive at a suitable generalization of the familiar fluctuation 
relations: 



compatible with the existence of macrostates with anomalous values in response 
functions, such as negative heat capacities Cy < 0, isothermal compressibilities K? < 
or isothermal magnetic susceptibilities xt < in those magnetic systems where this 
response function is expected to be nonnegative. 




(3) 




(4) 
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2. The proposal 

2.1. Notation conventions 

From the standard perspective of statistical mechanics, a system-surroundings 
equilibrium situation with several control parameters is customarily described by using 
the Boltzmann-Gibbs distributions |24j: 

dp BG (U,X\(3,Y) = z p — exp [-0 (U + YX)} Q (U, X) dUdX. (5) 

The quantities X = (V, M, P, Ni, . . .) represent other macroscopic observables acting 
in a given application (generalized displacements) such as the volume V, the 
magnetization M and polarization P, the number of chemical species Ni, etc.; with Y = 
(p, —H, —E, —fii, . . .) being the corresponding conjugated thermodynamic parameters 
(generalized forces): the external pressure p, magnetic H and electric E fields, the 
chemical potentials /ij, etc. In the sake of simplicity, we shall consider the Boltzmann's 
constant ks = 1 along this work, so that, in Eq.(j5]) represents the inverse temperature, 
= l/T. 

Conventionally, the system energy U is the most important physical observable in 
thermodynamics and statistical mechanics, and hence, it is always distinguished from 
the other macroscopic observables X. However, we shall not start from this distinction in 
the present approach, a consideration that allows us to deal with a more simple notation 
in our analysis. Let us consider the following convention for the physical observables 
(U,X) -> / = (J 1 ,/ 2 ,...) and (0,Y) ->• /? = (fit, fa,...) for the thermodynamic 
parameters, in a way that the total differential of entropy S can be rewritten as follows: 
dS = /3 (dU + YdX) dS = fadl 1 + (3 2 dl 2 + ... = frdP. We will use hereafter the 
Einstein's summation convention, which allows to rewrite the probabilistic weight of the 
Boltzmann-Gibbs distribution (JH]) as: 

u BG {I\0) = ^- ) e W {-0 l P). (6) 

2.2. Starting considerations on the distribution function 

At first glance, all that it is necessary to demand in order to extend fluctuation relations 
01]) is to start from a probabilistic distribution more general than the Boltzmann-Gibbs 
distributions (jSJ). From a mathematical point of view, a distribution that seems to 
describe the thermodynamical behavior of a system under a general influence of an 
environment is the ansatz: 

dp (I) = u(I)n (I) dl, (7) 

whose generic probabilistic weight oj (I) admits the Boltzmann-Gibbs form ()6| as a 
particular case. We shall shown in the next subsections that the hypothesis (jTJ), although 
simple, is sufficiently general to achieve our purposes. In fact, this expression is just a 
direct extension of the ansatz used in our previous work [10J: 

dp (U) =u(U)Q (U) dU (8) 
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to a more general case with several macroscopic observables. Let us now discuss under 
what general background conditions the mathematical form (j7]) can be supported by 
physical arguments. 

A particular and conventional situation is that where the system under study and its 
environment can be considered as two separable and independent finite subsystems that 
constitute a closed system in thermodynamic equilibrium. By admitting only additive 
physical observables (e.g.: energy, volume, particles number or electric charge) obeying 
the constraint It = I + I a, a simple ansatz for the distribution function is given by: 



where Q (I) and Q e (I e ) are the densities of states of the system and the environment, 
respectively, and W {It), the partition function that ensures the normalization 
condition: 



Here, T represents all admissible values of the system observables /. 

Although Eq.([9]) provides a simple interpretation for Eq.(J7J), it clearly dismisses 
some important practical situations where the additive constraint It = I + I e cannot be 
ensured for all physical observables involved in the system-environment thermodynamic 
interaction. Significant examples of unconstrained observables are the magnetization M 
or the electric polarization P associated with magnetic and electric systems respectively. 
A formal way to overcome such a difficulty is carried out by representing the density 
of state of the environment Q e in terms of the observables of the system I under study 
and the set of control parameters a of the given physical situation, Q e = Q e (I\a). 
Such an explicit dependence of the density of states of the environment on the internal 
state of the system follows from their mutual interaction, which leads to the existence 
of correlative effects between these systems. Thus, one can express the corresponding 
distribution function as follows: 



The subset r a of all admissible values of the physical observables I now depend on the 
control parameters a. 

Rigorously speaking, both separability and additivity represent suitable idealiza- 
tions of the real physical conditions associated with the large thermodynamic systems 
with a microscopic dynamics driven by short-range interactions. These conditions can- 
not be naively extended to systems with long-range interactions, which constitute a 
relevant framework for the existence of macrostates with anomalous response functions. 
As already mentioned in the introductory section, remarkable examples are the astro- 
physical systems, whose thermodynamic description is hallmarked by the existence of 
negative heat capacities C < 0. Although one cannot trivially divide a long-range inter- 
acting system into independent subsystems, an effective separability can be arisen as a 
direct consequence of the underlying microscopic dynamics in some practical situations. 




(9) 




(10) 




(11) 
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For example, it is physically admissible to speak about a separability between a globu- 
lar cluster and its nearby galaxy despite of the long-range character of the gravitational 
interaction. Since separability is, at least, an important condition to attribute a reason- 
able physical meaning to the description of a system as an individual entity, we shall 
demand this property as necessary requirement for the applicability of the ansatz (J7J). 

The existence of long-range interactions also implies the incidence of long-range 
correlations. This last feature, by itself, modifies in a significant way the thermodynamic 
behavior and the macroscopic dynamics of a physical system. If one considers the 
existing analogy with conventional systems at the critical points (where correlation 
length £ diverges, £ — > oo), the dynamics of long-range interacting systems should 
undergo the existence of some collective phenomena with a very slow relaxation towards 
the final equilibrium [9], as example, the so-called violent relaxation that appears 
during the collisionless dynamics of the astrophysical systems [33]. Under this last 
circumstance, it is usual that a long-range interacting system cannot be found in a final 
thermodynamic equilibrium, but in a long-living quasi- stationary state [251 EE]. The 
above dynamical features establish certain analogies among the long-range interacting 
systems with other physical systems with a complex microscopic dynamics such as 
turbulent fluids [34] and glassy systems [35] . 




According to the above reasonings, one cannot expect that the probabilistic weight 
uj(I) in Eq.((7]) can be always interpreted in terms of the density of states of the 
environment fl e as the cases associated with Eqs.flH]) or fTTTj) . since they demand the 
existence of a final thermodynamic equilibrium and the non-incidence of long-range 
interactions. Certainly, few systems in Nature are in absolute and final equilibrium, 
since it presupposes that all radiative materials would decayed completely and the 
nuclear reactions would have transmuted all nuclei to the most stable of isotopes. Thus, 
the hypothesis about the existence of a metastable equilibrium could be sufficient for 
the applicability of thermodynamics in many practical situations. Hereafter, we shall 
assume that the probabilistic weight uj(I) accounts for a general thermodynamic influence 
of an environment, which is found, at least, in conditions of meta- stability. This type 
of ad hoc hypothesis, that is, the using of generalized statistical ensemble with the form 
([7]), is usual in the statistical mechanics literature in the last decades. This kind of 
methodology has proved to be an useful alternative to describe thermodynamic features 
of systems with a complex microscopic dynamics. The reader can find an unifying 
framework of all these generalized distribution functions in the so-called Superstatistics, 
a theory recently proposed by C. Beck and E.G.D Cohen [36] . 

The specific mathematical form of the probabilistic weight depends, in general, 
on the internal structure of the environment, the character of its own equilibrium 
conditions, and the nature of forces driving the underlying system-environment 
thermodynamic interaction. Fortunately, the exact mathematical form of the function 
u (I) is unimportant in our subsequent development. Despite of its generality, we shall 
show that the internal thermodynamic state of the environment is fully characterized 
by a finite set of control thermodynamic parameters {3 U = (j3f, p% > ■ ■ ■) that constitute a 
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natural extension of the constant parameters ft = (Pi, P2, . . .) of the Boltzmann-Gibbs 
distributions (jSJ). Such a reduced description is analogous to the equilibrium situation 
described within the Gibbs' canonical ensemble (j2J), where the external influence is 
determined by the thermostat temperature T regardless its internal structure and 
composition, as long as the thermostat exhibited a very large heat capacity and 
guaranteed a thermal contact with the system under analysis. 



2. 3. Thermodynamic control parameters of the environment 

The direct way to arrive at a suitable extension of the thermodynamic control parameters 
P = (Pi, 02, . . .) of the Boltzmann-Gibbs distributions flBj) is achieved by appealing to the 
conditions of thermodynamic equilibrium such as the condition of thermal equilibrium, 
rpu _ rps^ foe condition of mechanical equilibrium, p^ = p s , the condition of chemical 
equilibrium, fi w = fi s , etc. As usual, the conditions of thermodynamic equilibrium can 
be derived as the stationary conditions associated with the most likely macrostate I of 
the distribution function p(I) = u (I) Q (I): 

sup {u (/) n (/)}=► A log {u (I) n(i)} = o, (12) 

which can be conveniently rewritten as follows: 

# (I) = PI (I) • (13) 
Here, /3£ are the system (microcanonical) thermodynamic parameters: 

derived from the coarsed grained entropy S (I) = log W (I) with W (I) — Q (I) 5c (5c 
is a small constant which makes the volume W dimensionless) , while the quantities 0f 
defined by: 

are the control thermodynamic parameters of the environment. This last definition 
constitutes a fundamental consideration in the present development. 

As expected, the environment control thermodynamic parameters (fT5l) drop to the 
constant parameters of the Boltzmann-Gibbs distributions (J6j) by replacing the 
general probabilistic weight u> (I) by the one associated with this last ensemble ubg (-0- 
Similarly, they drop to its microcanonical definitions (fT4l) when the environment is 
just a short-range interacting system in a final thermodynamic equilibrium, that is, 
when the probabilistic weight u is proportional to the states density Q e . However, 
these thermodynamic parameters turn effective in the other cases, overall, in those 
practical situations where the probabilistic weight u (I) corresponds to an environment 
in a metastable or quasi- stationary equilibrium. The possibility to extend this type 
of thermodynamic concepts to situations of metastability is not strange, since there 
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exist robust evidences supporting the introduction of concept of effective temperature 
in non-equilibrium systems (see in ref. [35J and references therein). 

A remarkable feature of the environment thermodynamic parameters ( fl5l) is their 
implicit dependence on the instantaneous macrostate I of the system, since these 
thermodynamic quantities are only constant for the particular case of the Boltzmann- 
Gibbs distributions. Consequently, it is almost a rule rather than an exception that 
the underlying system- environment interaction provokes the existence of non-vanishing 
correlated fluctuations among the environment thermodynamic parameters and system 
observables, (5j3^5P) 7^ 0. Such a realistic behavior of the interacting macroscopic 
systems is systematically disregarded by the use of Boltzmann-Gibbs distributions ([6j). 

The quantities act as control parameters during the system-environment 

thermodynamic interaction. A simple way to illustrate this role is by using the 
probabilistic weight u> (I) in the Metropolis Monte Carlo simulation [37] of the system 
under the present equilibrium conditions. The acceptance probability for a Metropolis 
move p ( I\ I + AJ) is given by: 



The fluctuations of the macroscopic observables satisfy the condition |AJ| <C |/| 
whenever the size N of the system is sufficiently large, allowing in this way to introduce 
the following approximation: 



and rewrite the acceptance probability (fl6l) in an analogous fashion to the usual 
Metropolis algorithm based on the Boltzmann-Gibbs weight ([6]): 



This last expression shows that the environment thermodynamic parameters ffl5l) 
provide a simple extension of control parameters of Boltzmann-Gibbs distribution (|6j). 
Latter, we shall employed this acceptance probability to perform numerical simulations. 

2.4- Derivation of some general fluctuation relations 

As already commented, the set T represents all admissible values of the system 
observables / = (J 2 , J 2 , . . . I n ). Let us denote by diA(I) = dA(I) /dP the first 
partial derivatives of a given real function A (I) defined on T. In general, the 
distribution function p(I) = of Eq.(J7J) should obey some mathematical 

properties associated with its physical significance. The are the following: 

CI. Existence: The distribution function p is a nonnegative, bounded, continuous and 

differentiable function on F. 
C2. Normalization: The distribution function p obeys the normalization condition: 




(16) 




(17) 






(19) 
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C3. Boundary conditions: The distribution function p vanishes with its first partial 
derivatives {dip} on the boundary dT of the set T. Moreover, the distribution 
function p satisfies for a < n + 1 the condition: 

lim \I\ a p(I) = 0, (20) 

|/|->oo 

when the boundary dT contains the infinite point {oo} of R n . 

C4. Existence of the second partial derivatives: The distribution function p admits 
second partial derivatives {didjp} almost everywhere on T, with the exception of a 
certain subset E cT with vanishing measure p (E): 

H(V)= [ p(l)dl = 0. (21) 

The mathematical conditions (C1-C3) constitute a direct extension of conditions 
(C1-C3) of our previous work to a case with several system observables (see in subsection 
3.3 of Ref.fTU]). Now, we have modified condition (C3) with the inclusion of the property 
( 1201 . as well as we have considered a fourth requirement, the condition (C4). In general, 
conditions (C3) and (C4) ensure the existence of the expectation values (P) and some 
relevant correlations functions such as (8P5P). 

Let us now introduce the quantity rji = r)i (J) = (3f (I) — (3- (I) defined as 
the difference between the i-ih thermodynamic parameters of the environment and 
the system, respectively. By taking into consideration definitions ( |T4l) and ( |T5i) . the 
thermodynamic quantity rji can be associated with a differential operator fji defined by: 

fji = -di (22) 

due to the validity of the mathematical identity: 

fj iP (I) = Vi p(I). (23) 

Its demonstration reads as follows: the partial first derivative of p can be rewritten as 
— dip = (—di log p) p and — di log p = —di log u — diS = fif — 0- = rji, with S = log W 
being the coarsed grained entropy with W = Q5c. 

The set of operators {fji} satisfy the following identities: 

[A{I),rj i ] = d i A{I), [7)i%] = 0, (24) 



with 



A, B 



AB — BA being the so-called commutator between two operators A and 



B. A remarkable particular case of the first commutator is obtained when A = P = P: 



I j ,Vi 



8j. (25) 



The consecutive application of two operators fji on the distribution function p leads to 
the following result: 

ViVjP = Vi (Vjp) = Vi (Vjp) = (~diVj + ViVj) P- (26) 
Notice here that condition (C4) is crucial for the existence of the partial derivatives d^j. 
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Now, let us combine the general mathematical properties (C1-C4) of the distri- 
bution function p (I) with the operational representation fl23|) of the thermodynamic 
quantity T]i to demonstrate the validity of the following expectation values: 

(t*) = 0, {mP) = H, (-fyfc + Wi> = 0, (27) 
which arise as particular cases of the following identity: 

(diA) = (Arji) (28) 

when A — 1, A — P and A = rjj, respectively. Here, A(I) is a differentiate function 
that satisfies the condition A(I) < \\I\ 2 when |/| — > oo if the set T is not finite. 

The demonstration start from the consideration of a set of arbitrary constants {a 1 }, 
which is employed to rephrase the integral: 



a 1 (Am) = J Aa\pdl (29) 
by using a vector J with components J 1 = a % p as follows: 

J Aa\pdl = J Aa%pdl = - J AV ■ Jdl. (30) 
By using the identity: 

AV ■ J = -V • (Afj + J ■ VA, (31) 
the integral (130]) can be rewritten as follows: 



AV ■ Jdl = - AJ ■ da + / J ■ VAdl. (32) 
Jdr Jr 

The surface integral over the boundary <9r vanishes as a consequence of condition (C3). 
This conclusion is evident when the set T is finite since J\gr = 0. Conversely, this 
vanishing takes place as a consequence of the limit: 

lim lirWW)! < lim A|/r +1 |p(/)l =0. (33) 

Identity fl28|) is obtained by considering a- 7 = for j ^ i and a 1 = 1 in the final 
expression: 

a? (Arij) = J J- VAdl = a j (djA) . (34) 

By taking into account the vanishing of (rji) and well as the definition of the 
correlation function (5r]i5P) = (rjiP) — (r/i) (P), the second relation of Eq.( )27|) can 
be conveniently rewritten as follows: 

(5 Vl 5P) = 51 (35) 

which is hereafter referred to as the fundamental fluctuation relation. By using this 
same procedure, the third relation of Eq.(l271) can be rephrased as follows: 

(diijj) = {SrjiSrjj) , (36) 
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which is hereafter referred to as the associated fluctuation relation. Let us now analyze 
the physical consequences of these mathematical results. 

The vanishing of expectation values (rji) in Eq.( l271) drops to the following 
thermodynamic relations: 

(P?) = (#> , (37) 

whose mathematical form explains by itself that the conditions of thermodynamic 
equilibrium are not only supported by the argument of the most likely macrostate 
(see Eq.( lT5|) in the previous subsection), but it can be rigorous expressed in terms 
of statistical expectation values. In fact, the present version ( 157)) of the thermodynamic 
equilibrium conditions has a more general applicability than the conventional version 
( 1T5j) . since the interpretation of the state of equilibrium through the most likely 
macrostate I is physically relevant for the large thermodynamic systems where thermal 
fluctuations can be disregarded. 

Conventionally, pairs of thermodynamic quantities such as energy U and 
temperature T, volume V and pressure p, magnetization M and magnetic field H, etc., 
are regarded as conjugated thermodynamic quantities. In this framework, commutative 
relation ( 1251) allows to refer to the pairs {P,f]i) as thermodynamic complementary 
quantities. Examples of thermodynamic complementary quantities are the energy U 
and the inverse temperature difference T]u = 1/T W — 1/T S between the environment 
(u) and system (s), or the generalized displacement X and the quantity r/x = 
/T w — Y s /T s defined in terms of the generalized forces (F s ,y w ) of the system and 
the environment respectively. The physical relevance of this new conceptualization 
relies on the fact that the thermal fluctuations of thermodynamic complementary 
quantities cannot be simultaneously reduced to zero. In order to shown this claim, 
let us first notice by considering Eq. (155l) that only the fluctuations corresponding 
to thermodynamic complementary quantities (P, r)i) are statistically dependent. By 
rewriting the fluctuation relations of complementary quantities through the known 
Schwarz inequality: 

(5A 2 ) (5B 2 ) > (5A5B) 2 , (38) 

one obtains the following remarkable relations: 

A Vt AP > 1, (39) 



where Ax = a/ (d~x 2 ) denotes the thermal dispersion of the quantity x. These inequalities 
obey the form of uncertainty relations rather analogous to the ones associated with 
quantum theories. In fact, Eq. (159l) is a direct extension of the energy-temperature 
uncertainty relation obtained in our previous work |10j : 

Af/A(l/T w - 1/T S ) > 1 (40) 

to other conjugated thermodynamic quantities such as the generalized displacement X 
and the generalized force Y: 

AlA(r/r - Y s /T s ) > 1. (41) 
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By generalizing our previous result, the physical significance of thermodynamic 
uncertainty relation fl39|) is the following: 

Claim 1: It is impossible to perform a simultaneous exact determination of the 
system conjugated thermodynamic parameter P and (3- = diS by using an experimental 



procedure based on the thermal equilibrium condition ( fff7| ) with a second system regarded 
here as a "measuring apparatus" : any attempt to reduce to zero the thermal uncertainties 
of the differences between their thermodynamic parameters Ar/j — > involves a strong 
perturbation on the conjugated macroscopic observable AP — > oo, and hence, P 
becomes indeterminable; any attempt to reduce this perturbation to zero, AP — > 0, 
makes impossible to determine the conjugated thermodynamic parameter [if by using 
the condition of thermal equilibrium (37\ ) since Arji — > oo. 

The present result constitutes a clear evidence about the existence of complementary 
relations in any physical theory with a statistical formulation, an idea postulated by 
Bohr in the early days of quantum mechanics [131 EE] with a long history in literature 
[T5l fl6l fT71 fT8l fT9l 1201 l2Tj . Although our results are intimately related with other 
approaches based on fluctuation theory [T5], [T71 120] . there exist subtle but fundamental 
differences distinguishing the present proposal from other attempts carried out in the 
past. The interested reader can see a more complete discussion about this subject in 
Ref.[38]. 

Let us now discuss some consequences of the associated fluctuation relation ( 1361) . 
The quantities 

Cij = (diVj) = ~ J P (I) didj log p (I) dl (42) 

are just the expectation values of the Hessian of the logarithm of the distribution 
function (I): 

P ij (I) = -d i d j logp{I), (43) 

which can be referred to as the probabilistic Hessian. Due to the positive definite 
character of the correlation matrix O = [O^-] with components = (5^57] j): 

aVOy = ((a%) 2 )>0, (44) 



as well as the applicability of the associated fluctuation relation fl36|) . Qj = Oij, the 
quantities Qj constitute the components of a positive definite symmetric matrix ( = [Qj] . 

The positive definite character of matrix [Qj] does not imply a positive definite 
character of the probabilistic Hessian ( )43l) for any macrostate / £ T. By admitting 
the existence of a certain region A w C V where this last condition does not hold 
for any admissible value of the internal control parameters that drive the form of 
the probabilistic weight the positive definite character of Qj indicates that 

macrostates belonging to the subset have always a poor contribution to the integral 
of Eq. (]42l . If the system size N is sufficiently large, the thermal fluctuations SP and 



| For example, the thermodynamic parameters (/3, Y) of the Boltzmann-Gibbs distribution ([5]), the 
total additive quantities It of distribution (j9|), or the formal parameters a of distribution 
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Srji turn negligible and the distribution function p (I) significantly differs from zero in a 
very small neighborhood around its local maxima. In this limit, all system macrostates 
I G A u turns inaccessible or thermodynamically unstable. Since the probabilistic Hessian 
(I43p depends on the particular mathematical form of the probabilistic weight u (I) 
that accounts for the external influence of the environment, it is possible to claim the 
following: 

Claim 2: The thermodynamic stability conditions and the fluctuating behavior of a 
given open system in (a quasi-stationary) equilibrium depend on the nature of external 
conditions which have been imposed to. 

This last conclusion appears to be self-evident and natural, but it contracts with 
some familiar equilibrium fluctuation relations such as the ones shown in Eq.(j3]), which 
seem to suggest (falsely) that system stability and fluctuating behavior only depend on 
its intrinsic properties: the response functions. However, such a conclusion is strictly 
correct for an equilibrium situation described by the Boltzmann-Gibbs distributions 
OH]). An illustrative example is the framework of application of the canonical fluctuation 
relation C = f3 2 (8U 2 ). The internal energy U fluctuates when the system is put in 
a thermal contact with a bath, but it remains constant when the system is put in an 
energetic isolation, and hence, there is no connection between the energy fluctuations 
and the heat capacity for the second situation. 

2.5. Gaussian approximation 

Fluctuation relations ( 135]) and f )36|) constitute two general mathematical results. The 
price to pay in order to arrive at a suitable extension of the familiar equilibrium 
fluctuation relations (jlj) is the consideration of a Gaussian approximation for thermal 
fluctuations, which appears as a licit treatment when the size N of the system under 
consideration is sufficiently large. 

Let us admit that the distribution function ([7]) only exhibits one (Gaussian) peak 
around the most likely macrostate J: 



with 5P = P — P. Since the thermal fluctuations 51 are small, the fluctuations 
of any continuous and different iable function A (I) can be expressed in a first-order 
approximation as follows: 



Formally, the thermodynamic parameters of the system /3? = diS (I) depend on the 
macroscopic observables /. Thus, their fluctuations S/3- can be expressed by using the 
fist-order approximation (f46l) : 




(45) 



5 A (I) ~diA (I) SP. 



(46) 



5Pt = d k Pt (I) SI k = H lk (I) 5I k , 



(47) 
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with Hij (7) = didjS (I) being the value of the entropy Hessian at the most likely 
macrostate. The substitution of Eq. (j47p into the fundamental fluctuation relation fl35l) 
allows to obtain the following expression: 

(5{3?5P) = 5\ + H ik (7) (5I h 5P) . (48) 

By using these same reasonings, the associated fluctuation relation ( 1361) can be rewritten 
as follows: 

d iVj (7) = (HH) • (49) 
Result P9l) evidences that the regions of thermodynamic stability are those 
macrostates where the probabilistic Hessian: 

i* ■ (7) = d iVj (I) = -d t d, {log u(l)+S (7) } , (50) 

is a positive definite matrix. Besides, Eqs. (1480 and ( ]4*9l) indicate that the system 
fluctuating behavior and its stability conditions depend on the external conditions 
which have been imposed to (Claim 2), which are characterized here by the correlations 
{5f3f5P) and the Hessian of the logarithm of probabilistic weight with opposite sign 
= —didjlogu (J). 

It is not difficult to verify that these identities provide the same qualitative 
information about the system thermodynamic stability. For example, the equilibrium 
conditions associated with the Boltzmann - Gibbs distributions (EJ) follows from the 
imposition of constraint 8/3£ = into Eq. (f48|) . Thus, the correlation matrix C 1 ^ = 
(5I k 5P') is given by the inverse matrix of the entropy Hessian (I) with opposite 
sign: 

-H ij (I)=C ij , (51) 

which is a well-known result of classical fluctuation theory |24j. Since the correlation 
function C 1 ^ is always a positive definite matrix: 

a iaj C ij = /( ai 5P) 2 \ > 0, (52) 

the peaks of the Boltzmann- Gibbs distribution ([5]) will never be located in those regions 
where Hij (I) is a nonnegative definite matrix, that is, where the entropy S is not a 
locally concave function. Thus, these macrostates turn inaccessible in the limit of large 
system size N, so that, all them constitute the unstable region = A BG of this 
particular ensemble. Clearly, this is the same conclusion derived from Eqs.( H9|) and 
(|50|) . As a by-product of the present reasonings, we have arrive at the following: 

Claim 3: The direct observation of unstable macrostates I G A BG of the 
Boltzmann- Gibbs distributions (TJJ) is only possible when (5(3^5P) ^ 0, that is, by 
considering an equilibrium situation where the internal thermodynamic state of the 
environment is affected by thermodynamic interaction with the system. 

There exist within the Gaussian approximation a fluctuation relation which has no 
counterpart within the Boltzmann-Gibbs distributions since it involves the correlation 
matrix = (5(3^5(3^. By using the first-order approximation: 

5ft = djff (I) 5P = Fij (I) 5P (53) 
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(54) 

Ha (I). By 



In the sake of simplicity, let us adopt the notation Fy = Fy [Tj and Hi 
applying two times the fluctuation relation (HHl) and the approximation 



(5p?6P) = F im {5I m 6P), 



(55) 



one obtains: 



(8/3?8/3?) = Fij + E i3 + H im C mn H, 



(56) 



This result can be conveniently rewritten by performing the matrix product with the 
correlation matrix G %3 = (8P8P): 



It is possible to recognize that the right hand side of the above equation is simply the 



term M k M k , where M{ =51 + Hn-C^ denotes the correlation matrix Mf = (8/3" 6 ft) 
of Eq. pSj) . Thus, we have derived a complementary fluctuation relation: 



very difficult to measure in a given experimental situation, since it demands to 
perform simultaneous measurements of the system observables I and the thermodynamic 
parameters of the environment (3 e . The complementary fluctuation relation ( 1591) allows 
to obtain in an indirect way these fluctuating quantities by performing independent 
measurements of the fluctuations of the system observables I and the environment 
thermodynamic parameters (3 e . 

2.6. Extended equilibrium fluctuation- dissipation theorem 

For convenience, let us re-derive here the fluctuation relations associated with the 
Boltzmann-Gibbs distribution ([6]). The direct way to carry out this aim is by starting 
from the partition function Z (/?): 



Let us denote by d l A(/3) = dA(/3) /d(3i the first partial derivatives of an arbitrary 
function A (/3) defined on the thermodynamic parameters {3. It is easy to verify the 
validity of the following identities: 



r> /~tkj 17 s-ikj i tt /~ikj i tj r~imn tt s~<kj 

= 8j + 2H ik C kj + H m C mn H nk C k K 



(57) 





(60) 




(61) 
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The substitution of the first relation into the second one yields: 

X% G = C ii = {5I i 8P), (62) 
where the symmetric matrix Xbg- 

X%o - ~ d -§± (63) 

will be referred to as the canonical response matrix. Since the correlation matrix 
C" J = (8P8P) is a positive definite matrix, the canonical response matrix Xbg ls a ^ so a 
positive definite matrix. 

The identity ( )62l) is a rigorous result within the Boltzmann-Gibbs distributions fl6]). 
Notice that an analogous expression can be derived within the Gaussian approximation 
from the fluctuation relation ( )48|) after imposing the restriction 8(3f = 0, that is, Eq.tlBTT). 
which can be rephrased as: 

= (8P8P) , (64) 

with x l J being the inverse matrix of the entropy Hessian with opposite sign evaluated 
at the most likely macrostate: 

X* = -U" (I) . (65) 

Comparison between Eqs.( l62l and f l64"l) leads to establish the relationship Xbg ~ X % s > 
which allows to refer to x l J as the microcanonical response matrix of the system. 

The microcanonical response matrix x % i defined above is a purely microcanonical 
quantity, that is, it is an intrinsic thermodynamic quantity that can be attributed to an 
isolated system. Actually, the physical meaning of such a response matrix is unclear. In 
general, a response function characterizes the sensibility of the system under the external 
control, and consequently, such a concept is only relevant for an open system. While such 
a physical relevance is satisfied by the canonical response matrix Xbg^ applicability 
is only restricted to the case when the thermodynamic influence of the environment 
corresponds to the Boltzmann-Gibbs distributions ([6]), where the thermodynamic 
parameters are constant parameters instead of fluctuating quantities. 

According to these last reasonings, the concept of response matrix x % ^ should be 
extended in order to be also applicable in a more general system-environment equilibrium 
situation. Since the physical observables / and the (control) thermodynamic variables 
/3 U undergoing, in general, thermal fluctuations, a simple operational definition for the 
concept of response matrix Xu could be introduced in terms of expectation values as 
follows: 

d(P) d(P) 

Here, we have took into account the equilibrium conditions (fif) = (j3f). 

Clearly, definition (16 6 j) contains the canonical response matrix Xbg as a particular 
case when 8^ — =r> fif — ~ const. The response matrix x% defined in this way is 
more realistic and general than the one associated with the canonical response matrix 
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( 16"5j) . The canonical response matrix Xbg presupposes that the control parameters 
{Pi} have no fluctuations. Such an exigency is clearly an idealization since any real 
measuring process always involves a perturbation of the thermodynamic state of the 
measuring apparatus. For example, both the temperature indicated by a thermometer 
and the pressure indicated by a barometer are affected by the thermodynamic influence 
of the system under study. In practice, the true values of temperature T and pressure p 
employed to calculate the heat capacity C p at constant pressure and the isothermal 
compressibility Kt of a given system are actually expectation values, (T) and (p). 
Therefore, the practical operational definition of these concepts reads as follows: 



\d(T)J {p} (V)\d(p) /m 
where we have taken into consideration the stochastic character of the energy U, the 
volume V and the entropy S = S (U,V). 

As already discussed, the expectation values (P) and (/??) are given within the 
Gaussian approximation by the corresponding values at the most likely macrostate I, 
(P) = P and ((3?) = /3- = /3- (7), which allows to rewrite Eq. fl66|) as follows: 

BP 

Xt = -9fi = ft ■ (68) 

Claim 4: The response matrix x l J °f (I6"6"j) could be identified with the 
microcanonical response matrix x l J f|65|) at the most likely macrostate I as long as 
the Gaussian approximation applied. 

Equilibrium fluctuation relation (I64p is inapplicable within the anomalous region 
A BG of the Boltzmann-Gibbs distributions because of the microcanonical response 
matrix x l J turns a non positive defined matrix. 

Claim 5: The unstable region of the Boltzmann-Gibbs distributions A 
corresponds to macrostates with anomalous values in the response functions, such as 
macrostates with negative heat capacities at constant magnetic C, p < or negative 
isothermal compressibilities Kt < in a fluid system, or macrostates with negative 
heat capacities at constant magnetic field Ch < or isothermal magnetic susceptibilities 
Xt < in a paramagnetic system. 

As already shown, the presence of macrostates with anomalous values of response 
functions is incompatible with the imposition of the restriction 5(5" = into the 
fluctuation relation ( )48l) . However, this same fluctuation relation can be rewritten to 
obtain an extension of Eq. flM]) . which is expressed here in the following: 

Theorem: The system fluctuating behavior C 1 ^ = ( x 5I k 5P^ is determined within 
the Gaussian approximation by its response x % i — ~ d (P) /d((5f) and the incidence of 
correlative effects Mf = (5(3%5P) that appear as a consequence of the external influence 
of the environment through the following fluctuation relations: 

xii = C«+xi k Ml. (69) 

Eq. (l69l) constitutes the central result of this work. Since the fluctuation relation 
(1621) is simply the equilibrium version of the so-called fluctuation- dissipation theorem, the 
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expression (IBTJj) can be unambiguously referred to as the extended equilibrium fluctuation- 
dissipation theorem. The same one accounts for the system fluctuating behavior 
beyond the equilibrium thermodynamic conditions associated with the Boltzmann-Gibbs 
distributions (EJ in a way that it is compatible with the existence of anomalous response 
functions. As already discussed, its validity relies on the applicability of the Gaussian 
approximation, which arises as a licit treatment in the limit of large system size N. 
As expected, the present theorem comprises all those preliminary conclusions obtained 
along this work, with the exception of claim 1, which was obtained from a more general 
equilibrium fluctuation relation, Eq.( l35|) . 

Let us finally obtain the form of the extended equilibrium fluctuation-dissipation 
theorem fl69|) in terms of the ordinary thermodynamic variables and potentials of 
statistical mechanics and thermodynamics: (U, (3) and (X, Y) and the Enthalpy G = 
U + YX. The known relations: 

-dp log Z (P, Y) = (G) , (3 (X) = -dy log Z (f3, Y) 
d 2 log Z (P, Y) = (6Q 2 ) , d 2 log Z (0, Y) = (3 2 (6X 2 ) , (70) 
df,dylogZ(P,Y) = p(5Q5X), 

allow to re-express Eq. (j62l) in a compact matrix form as follows: 

11 = 7, (71) 

where the response and the correlation matrixes TZ and J 7 are given by: 

d p {G) d p {p{X)) 

dy (G) (5dy (X) 

(6Q 2 ) f3 (6Q5X) 
(3 (6X5Q) (3 2 (5X 2 ) 

It is important to keep in mind that the Enthalpy fluctuation 5G within the Boltzmann- 
Gibbs distribution (j5J) is amount of heat 5Q exchanged between the system and its 
environment at the equilibrium, 5G = 5U + Y5X = 5Q. Such a relationship does 
not hold in general, since SG = SQ + X5Y when the fluctuations 5Y are taken into 
consideration. The matrix form of the fluctuation-dissipation relations (162|) in terms of 
the conjugated thermodynamic variables (U, (3) and (X; £) with £ = (3Y reads as follows: 

X = C, (73) 
where the response and the correlation matrixes x an d C are given by: 
x= _( d P (U) d p (X)\ 

x { d,(u) o,(x) y [ ' > 

r _( (SU 2 ) (5U5X}\ 
U ~[(5X5U) (5X 2 ) J' 

One can verify the existence of the following transformation rules between these 
two representations: 

jr = TCT T , K = TxT 71 , (76) 



n 



(72) 
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where the transformation matrix T is given by: 

Finally, one can rewrite the compact matrix form of extended equilibrium fluctuation- 
dissipation theorem f )69|) in these equivalent representations as follows: 

X = C + X M^'ll = J r + TIM (78) 

where the correlation matrixes M and A4 are given by: 



M 



M 



(6^511) (5P U 6X) 
(8FSU) (SCSX) / 

(5P W 6Q) PQPPSX) N 



{5Y W 5Q) f3{8Y"8X) 

which are related by the transformation rule: 

M = {T T )~ l MT T . (80) 

As expected, the quantities (P u , Y u ) denote the thermodynamic parameters of the 
environment. 

3. Examples of application 

3.1. Application to fluid and magnetic systems 

Direct application of Eqs. (172p and (!79|) to the thermodynamic variables of a fluid system 
with Enthalpy G = U + pV leads to the following expressions of the matrixes 71, J 7 and 
M: 



K 



T 2 C P V {Ta p - 1) 
(dG/dp) T PVK T 

{5Q 2 } (6Q5V) 

(5V5Q) P 2 (5V 2 ) 

(5p w 5Q) p {5p w 5V) 



oTZL t;; i, (si) 



M 



where C p = (dG/dT) p and a p = V^ 1 (dV/dT) p are the heat capacity and the thermal 
expansibility at constant pressure, while K? = —V~ x (dV/dp) T is the isothermal 
compressibility. The symmetry of the response matrix in Eq. (18 1 j) leads to the 
thermodynamical identity: 

dG\ _fdV\ , . 
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The corresponding expressions for a magnetic system with Enthalpy G = U — HM 
are directly derived from the previous results by substituting (V,p) — » (M, —H): 



where Ch = {OG/dT) H is the heat capacity at constant magnetic field and xt — 
(dM/dH) T the isothermal magnetic susceptibility. The symmetry of the response 
matrix in Eq. fl83|) leads to the thermodynamical identity: 



3.2. Thermodynamic stability of macro states with anomalous response functions 

A direct consequence of extended equilibrium fluctuation-dissipation theorem (1751) are 
the conditions of thermodynamic stability, which can be compatible in this framework 
with the existence of anomalous response functions. In the present subsection, let 
us discuss two simple application examples involving macrostates with negative heat 
capacities and a numerical testing during the study of macrostates with anomalous 
values in the isothermal susceptibility of a paramagnetic system. 

3.2.1. Thirring's constraint: Let us first consider an equilibrium situation involving 
only the conjugated thermodynamic pair energy-temperature. Since 5X = =>■ 
5Q = SU, Eq. (1751) drops to the generalized equilibrium fluctuation-dissipation relation 
involving the heat capacity ([3]) already obtained in our previous work |10]. A simple 
way to implement the existence of non-vanishing correlated fluctuations (5/3 w SU) is by 
considering a thermal contact with a bath (B) having a finite heat capacity Cb > 0. 
Clearly, the bath experiences a temperature fluctuation 5T B when the system absorbs 
or releases an amount of energy 6U, 5Tb = —5U/Cb =>■ (Sf3sSU) ^ 0. Eq.(j3D can be 
rephrased as follows: 



after considering the Gaussian approximation 6(3b = (3 2 5U/Cb (here k B = 1). Since 
the expectation value of the square dispersion of the energy fluctuations (5U 2 ) is 
nonnegative, the thermodynamic stability of macrostates with negative heat capacities 
C < demands the applicability of the following condition: 




(83) 




(84) 




(85) 



CC B 



> o < C B < \C\. 



(86) 



C + C B 
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This result is same constraint derived by Thirring in the past As already 

commented, the existence of anomalous macrostates with negative heat capacities is 
a distinguish feature of many nonextensive systems. However, since the derivation of 
Thirring's constraint ( 186|) demands the existence of a thermal contact, such a stability 
condition could be only applicable to short-range interacting systems whose sizes are 
large enough in order to dismiss the surface energy involved in their thermodynamic 
interaction. 

3.2.2. Stability condition in an astrophysical situation: Let us now consider a 
paradigmatic situation with long-range interactions: the astrophysical system composed 
by a globular cluster (C) and its nearby galaxy (G). Let us denote by Cq and Cq their 
respective heat capacities, being \Cc\ <C \Cq\- At first glance, the galaxy could be 
taken as a bath with a very large heat capacity, which naively suggests that the thermo- 
statistical description of the globular cluster could be performed with the help of the 
Gibbs' canonical ensemble fl2]). In this context, however, the known exigence of a thermal 
contact fails as a consequence of the long-range nature of Gravitation. 

Astronomers known since the late 1800s that the attractive character and power- 
law 1/r 2 dependence of gravitational interaction favor a concentration of stars towards 
the inner regions of any astrophysical structure and the increasing of their kinetic 
energies while a decreasing of the total energy |40j . a tendency stimulating the 
development of distribution profiles with a core-halo structure and negative heat 
capacities. Consequently, it is a rule rather than an exception the fact that the globular 
cluster and its nearby galaxy simultaneously exhibit macrostates with negative heat 
capacities, Cc < and Cq < 0. Since Thirring's constraint ( 1861) cannot be satisfied, 
and overall, the same one is inapplicable in the present situation, the energy associated 
with the thermodynamic interaction between the globular cluster and the galaxy should 
play a fundamental role in the thermodynamic stability. 

To obtain a general idea about how works the above mechanism, let us denote by 
Wq-c the energy contribution of the total energy Ut = Ug + Uc + Wg-c accounting for 
the potential energy Vq-c associated with the gravitational interaction between these 
systems, as well as the kinetic energy contributions of the collective translational motion 



In general, Wq-c should be significantly smaller than the internal energy of the 
galaxy Uq, but it could be comparable to the internal energy of the globular cluster 
Uc, | Wg-c | ~ Wc\ <C Wg\- Strictly speaking, the total energy Ut changes as a 
consequence of gravitational influence of other neighboring astrophysical structures 
and the evaporation of stars. However, the composed system globular cluster-galaxy 
can arrive at a certain quasi-stationary evolution under its own relaxation mechanisms 
(collisions, violent relaxation, etc. [33]) and the incidence of stars evaporation, where 
Ut changes at a very slow rate, Ut — 0. In these conditions, it is possible to speak about 



Wg-c = Tg-c + Vg-c- 
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the existence of a certain metastable thermodynamic equilibrium between the globular 
cluster and its nearby galaxy, a fact that it is manifested in the observation of universal 
profiles such as the Vaucouleurs' law in the surface brightness of the elliptical galaxies 
and the Michie-King profiles for the distribution of stars in the globular clusters |JT]. 

Due to the coupling between collective and internal degrees of freedom, the energy 
contribution W G -c wm depends on the internal states of the globular cluster and 
its nearby galaxy, so that, it can be expressed by their respective internal energies 
(U G , Uc) as well as some other collective variables a (their total masses, relative 
separation, etc.), Wq-c — W(U G , Uc', cr). Consequently, the cluster and galaxy undergo, 
in general, a nonlinear energetic interchange during their thermodynamic interaction, 
AU G 7^ —AUc, since part of the energy transferred from the galaxy AUq towards the 
cluster is used to modify the interaction energy Wg-c- Such a nonlinear interchange 
contracts with the usual thermal contact between a Gibbs' bath (B) and an extensive 
system (S), where AUb = — AC/5. 

As a consequence of the separability as well as the small number of degrees of 
freedom of collective motions j\f^^ Uve [ n comparison with the internal degrees of 
freedom of the galaxy N™ ternal an d the globular cluster N™ ternal ; the distribution 
associated with the present situation can be approximately given as the product of 
their respective densities of states Qg an d &c'- 

p(Uc)^AQ g (Ug)^c(Uc), (88) 

where the internal energies Uq and Uc obey the constraint: 

U T = U G + U C + W G -c ~ const, (89) 

with A being a normalization constant. Since the probabilistic weight co(Uc) ~ & G (U G ), 
the environment inverse temperature (see definition in Eq. (fT5]) ) is given by: 

dS G (U G ) 



dU 



c 



= G a, (90) 

u T 

where the subindex Ut means that the partial derivation is taken at constant 
total energy. In this last expression, (3 G = dS G (U G )/dU G represents the ordinary 
microcanonical inverse temperature of the galaxy and a = a(Uc)' 



dU G 

a = — 



dU 



c 



^ l + dw G - c 



Wc 



= a(U c ), (91) 



a dimensionless factor that accounts for the underlying nonlinear energetic interchange. 
Clearly, the inverse temperature ( 190]) is an effective concept that characterizes how the 
globular cluster fills the thermodynamic influence of the galaxy. Since in general 
the condition of thermal equilibrium j3 w = ftc implies that the microcanonical inverse 
temperatures of the galaxy and the globular cluster are different, 0g 7^ As*- 

Since the heat capacity of galaxy C G is extremely large in comparison with the 
heat capacity of the globular cluster C G , the inverse temperature /3 G can be taken as a 
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constant parameter. Thus, we can limit to consider the thermal fluctuations of effective 
inverse temperature /3 W associated with the nonlinear energetic interchange: 

5/3" ~ f3 G a'(U c )6U c . (92) 
By substituting this last approximation into fluctuation-dissipation relation of Eq.(j3j) 
and considering the condition of thermal equilibrium (3 W = flc, one obtains: 

RS^ ^W- (93) 

Since (8Uq) > 0, one arrives at the following condition of thermodynamic stability: 

PcCc 1 + arV > 0. (94) 

According to the present result, the globular cluster with Cq < can be found in 
a metastable thermodynamic with its nearby galaxy provided that the underlying 
nonlinear energetic interchange associated with the interaction energy Wq-c ensured 
the applicability of stability condition (p2 



3.2.3. Anomalous response in paramagnetic systems: Besides the long-range 
interacting systems as the astrophysical example previously discussed, the existence of 
macrostates with negative heat capacities is a relevant feature indicating the occurrence 
of a discontinuous (first-order) phase transition in a mesoscopic short-range interacting 
system [8]. The observation of anomalies in discontinuous phase transitions is not 
restricted to the heat capacity, but other response functions can also exhibit anomalous 
values too [32J. In fact, their presence could be the origin of the sudden change 
experienced by a given physical observable X with a small varying of its conjugated 
control parameter Y below the critical temperature T c . 

A typical example is found in a magnetic system between the magnetization M 
and the external magnetic field H . The sudden inversion of magnetization M \- > —M 
with M/0, observed in turns to the value H = below the critical temperature T c , 
is a consequence of the existence of unstable diamagnetic states xt < 0. Such states 
are anomalous for magnetic materials undergoing a paramagnetic-ferromagnetic phase 
transition, in those equilibrium situations where the fluctuation-dissipation relation 
Xt = P {SM 2 ) is supposed to be relevant, and the inequality xt > arises as the 
ordinary stability condition. 

Let us consider again a magnetic system with internal energy U and magnetization 
M. As already shown in this work, the corresponding form of the generalized equilibrium 
fluctuation-dissipation relations 71 = F+TZM. is specified by matrixes of Eq. (l83]) . In the 
sake of simplicity, let us restrict to a particular equilibrium situation where temperature 
of environment /3 U shows a constant value 0, but the external magnetic field H u acting 
on the magnetic system undergoes non- vanishing thermal fluctuations (jSH^SM) ^ 0. 

The simplest way to implement the above behavior is by admitting the existence 
small thermodynamic fluctuation of the magnetic field if w in turn to a constant value 
H coupled to the total magnetization of the system: 

H w = H - X5M/N, (95) 
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where N is the system size and A a suitable coupling constant. The general equilibrium 
fluctuation-dissipation relation involving the isothermal magnetic susceptibility xt '■ 

P XT = p 2 (SM 2 ) + (T [dM/dT) H — M) (3 {5^5M) -(3 2 X t (SH^SM) (96) 

drops to: 

P XT = p 2 <5M 2 ) - (3 2 X t (5H"5M) (97) 
in the present equilibrium situation, which can be rewritten as follows: 

^ M2 > = TT^7F < 98) 

Since the standard deviation of magnetization is nonnegative, (5M 2 ) > 0, the anomalous 
macrostates with negative values of isothermal magnetic susceptibility xt < are 
thermodynamically stable when the condition: 

A + N/xt > (99) 

holds. Eq.(l98l) can also be used to obtain the dispersion of the external magnetic field 

< M ^ W >^ r ^ F (ioo) 

Since xt grows with iV in general as xt oc N, the dispersions of the magnetization 
and the external mag netic field behave as (5M 2 ) cx N and ((dH^) 2 ) oc l/N. When 
the system size N is large enough, the fluctuations of the external magnetic field 
are very small. From the viewpoint of thermodynamics (where thermal fluctuations are 
dismissed due to the large size of conventional thermodynamic systems), the present 
equilibrium does not differ from conventional situations where the external magnetic 
field is assumed to be constant. However, it is easy to see that both the fluctuating 
behavior accounted for expressions fl98|) and fllOOp . as well as stability condition (j99j) 
crucially depend on the external conditions which have been imposed to, which are 
specified here by the coupling constant A. 

According to the stability condition ( 199]) . macrostates with negative isothermal 
susceptibilities can be found in a stable equilibrium with an appropriate choosing of 
the coupling constant A. Even, fluctuation relation fl98|) could be also used to obtain 
the isothermal magnetic susceptibility per particle xt = Xt/N in terms of fluctuating 
behavior of magnetization: 

XT= ww) N= ^T' (101) 

where o~ 2 m = (5M 2 ) /N represents the thermal dispersion of magnetization. 

Certainly, we unknown how the external magnetic influence expressed in Eq. fl95l) 
could be implemented in a real experimental situation. Fortunately, their consequences 
are very easy to check with the help of Monte Carlo simulations. Let us consider the 
known 2D Ising model on the square lattice L x L with periodic boundary condition: 

U=-J2siSj, M = J2 s h ( 102 ) 

(V) 1 
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Figure 1. Isotherms of 2D Ising model with L = 25 derived from Metropolis Monte 
Carlo simulations by using the acceptance probability (|103[) : Panel a) magnetic field 
(H") versus magnetization (to); Panel b) inverse isothermal magnetic susceptibility 
X? 1 versus magnetization (to). The labels in both panels indicate the value of j3 for 
the corresponding isotherm. 



where the sum (ij) considers nearest-neighbor interactions only, and the spin variables 
Si = ±1. The above equilibrium conditions can be easily implemented by using a 
Metropolis method with acceptance probability p = p (U, M \ U + AU, M + AM): 

p = min {exp [-fiAU + 0H w AM] } , (103) 

which is a particular expression of Eq.(fl8l). By denoting m = M/N the magnetization 
per particle, the external magnetic field in this study is given by = H + A (m — m), 
where fh and H are suitable estimation of the expectation values (m) and (H^). Our 
interest is to obtain the isotherms of 2D Ising model within anomalous regions with 
Xt < 0. Thus, the values of parameters (ff, m) can iteratively provide by using the 
isothermal magnetic susceptibility per particle \t °f a previous simulation throughout 
the expression: 

H i+1 = Hi + (xt); 1 (m i+1 - ffn) , (104) 
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Metropolis Monte Carlo simulations: 
-»-with (&H°5M)*0 (the present method) 
-•-with <8W"SA/}=0 (BG) 




(m) 

Figure 2. Isotherms for L = 25 and f3 — 0.46 obtained from Metropolis Monte Carlo 
simulations based on the Boltzmann-Gibbs distributions (BG) where H u = H ~ 
const — > (SH^SM) = and the present Metropolis Monte Carlo simulations with 
(5H"5M) ^ 0. 



where the step Am = mj +1 — rrij should be small. The initial value m is assumed 
to be the average of magnetization obtained from an ordinary Metropolis Monte Carlo 
simulation with H 1 ^ = Hq ~ const (A = 0) far enough from the unstable region with 
Xt < 0, m = (m). 

Although any real value of coupling constant A satisfying stability condition ( 1991) is 
admissible, it is desirable to impose some criterium in order to reduce thermal dispersions 
of the system magnetization ( 1981) and the external magnetic field f llOOp . Our interest is 
to obtain the isotherms in the H — M diagram, and hence, an optimal value for A can 
be chosen in order to minimize the total dispersion a 2 = a 2 H + cr^, where: 

a 2 H = N{{5H"f) (105) 

is the thermal dispersion of the external magnetic field. The analysis yields: 

A op = Va 2 + 1 - a, (106) 

where a is the inverse of the isothermal magnetic susceptibility per particle, a = x^ 1 , 
whose value can be estimated from the previous simulation as the case of parameters 
(H,mJ in Eq. (I104j) . Notice that this optimal value guarantees the validity of stability 
condition (1991 . 

Results of extensive Monte Carlo simulations by using the procedure explained 
above are shown in FIGjlJ We have used a square lattice with L = 25 and considered 
n = 10 6 iterations for each calculated point. This results clearly revealed the presence 
of anomalous diamagnetic states xt < for inverse temperatures f3 above the critical 
point ftc =25 ~ 0.41. We only show here macrostates with (m) > due to the underlying 
symmetry M -> —M and H — >■ —H of the Enthalpy G = U - HM. 

As clearly shown in FIGJ2J while the 2D Ising model isotherms are able to describe 
a backbending in the H — M diagram under a magnetic influence obeying stability 
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condition f !99p . they merely undergo a sudden jump in magnetization in turns to the 
point H = under the influence of a constant magnetic field, that is, a situation with 
A = 0. As already commented, such a sudden jump is interpreted as the occurrence of a 
discontinuous phase transition in conventional statistical mechanics, whose existence is a 
direct manifestation of the presence of an unstable branch with negative susceptibilities 
Xt < 0. Our Monte Carlo simulations show that such an anomalous region can turn 
accessible with an appropriate selection of the external experimental conditions. The 
present results can be extended to the case of fluid systems. In fact, the behavior of the 
2D Ising model isotherms is fully analogous to the well-known backbending of Van der 
Waal gas isotherms in the p — V diagram, which is associated with macrostates with 
negative isothermal compressibilities < [42J. 

4. Final remarks 

We have obtained in this work a generalization of equilibrium fluctuation-dissipation 
relations compatible with the existence of macrostates with anomalous response 
functions, Eq. fl69|) . Such a development has been carried out by starting from a 
generic extension of the Boltzmann-Gibbs distribution and its control thermodynamic 
parameters (/3,Y) that characterize the external thermodynamic influence of an 
environment, Eqs.([7|) and (fl5l) . Since the environment thermodynamic control 
parameters turn, in general, fluctuating quantities, a new operational definition (1661) 
for the concept of response function has been also considered. Remarkably, the physics 
behind the generalized equilibrium fluctuation-dissipation relations ( )69|) is quite simple. 
The key ingredient relies on the that the environment thermodynamic parameters 
{Pi} undergo in general correlated fluctuations with the system observables {P} as 
a consequence of the underlying thermodynamic interaction. Such a realistic behavior, 
systematically disregarded by the using of Boltzmann-Gibbs distributions ([6]), leads 
to the consideration of a new correlation matrix M. in the equilibrium fluctuation- 
dissipation relations, Eq. (!79|) . which characterizes the existence of these correlative 
effects. The new theorem explicitly expresses that the system fluctuating behavior and 
its thermodynamic stability crucially depend on the nature of the external influence of 
the environment. 

As example of applications, we have obtained particular expressions of the 
generalized equilibrium fluctuation-dissipation relations for a fluid and magnetic 
systems. After, we have derived two criteria of thermodynamic stability accounting 
for the existence of negative heat capacities: Eq. flBBl . for systems with short-range 
interactions ensuring the existence of a thermal contact] and Eq.( !94l) . for a metastable 
equilibrium involving long-range interactions between a globular cluster and its nearby 
galaxy. Finally, we have arrived at a stability condition compatible with anomalous 
isothermal susceptibilities associated with the occurrence of a discontinuous phase 
transition in magnetic systems, Eq. fl99|) . whose consequences have been also numerically 
tested throughout Monte Carlo simulations. 
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We would like mention some open questions deserving a special attention in future 
works: 

• The usual equilibrium fluctuation-dissipation relations ( I7TT) constitute a particular 
expression of the so-called fluctuation- dissipation theorem: 



relating correlations and response functions in linear response theory [12]. The 
present generalization of equilibrium fluctuation-dissipation relations inspired a 
corresponding extension of this last fundamental result. 

• The numerical study of the 2D Ising model clearly shows that the present ideas are 
very useful in order to implement some new Monte Carlo algorithms inspired on 
statistical mechanics [12j in order to access to macrostates with anomalous response 
functions. 

• At first glance, rigorous fluctuations relations (|35l) and (|36l) clearly exhibit a 
tensorial character, which suggests a direct connection of the present approach with 
some geometric formulations of fluctuation theory developed in the past [221 123] . 
A preliminary analysis of such a question is presented in a recent paper [TTJ . 
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